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Ground-state energies are investigated in a many-quark model with pairing interactions, 
which has the su(4)-algebraic structure. Exact eigenstates in the boson realization method 
are constructed by imposing a color-singlet condition developed in the previous paper. An 
interaction term breaking the stt(4)-dynamical symmetry plays an important role to deter- 
mine the ground state. As a result, a quark-pairing state or a quark-triplet state as a nucleon 
is realized with a certain value of a variable which is regarded as an order parameter. In 
addition to the parameter regions in which these ground states are realized, it is shown that 
there are two transition regions between the quark-pairing and the quark-triplet states with 
different values of the order parameter. 

o 

§1. Introduction 

The Bonn quark model is an interesting model, which was first introduced with 

a purpose of describing the nucleon and the Z\-resonance as quark-triplet states.^ 1 

' The original Bonn quark model has the following Hamiltonian 

t-h ; 

H = - ^2( C *2m c 3m c 3m'C2m' + ct m C* lfji C im c 3m> + C*m c 2fh c 2m 'Clm') , (1-1) 
. | mm' 

where c* m and Cj m are quark creation and annihilation operators with color i and 
the angular momentum quantum number of the single quark level, m = —j 8 , —j s + 
1 , ••• ,j s . Here, c* iih = (— l)- ?s ~ m c*_ m . This Hamiltonian represents the quark- 
pairing interaction which, in general, leads to the color instability. Namely, it is 
possible that a color superconducting state is realized. It has been remarked that 
the original Bonn quark model has a dynamical szi(4)-symmetryJ^ Thus, a su(4)- 
symmetry breaking term is introduced in which the color su(3)-symmetry is retained. 
Namely, a su(4)-symmetry breaking interaction proportional to the s-u(3)-Casimir 

~2 

operator, Q , is introduced, which represents a particle-hole-type interaction in terms 
of the quark shell model: 

H m = H + X Q 2 ■ (1-2) 

This model is called the modified Bonn quark modelP 1 In the previous paper, in 
RefJ2D, which is hereafter referred to as (A), exact eigenstates were investigated by 



typeset using -pTpf^ ds (Ver.0.9) 
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the method of the boson realization. In Ref.3), which is hereafter referred to as 
(B), the exact eigenstates with single-quark, quark-pair and quark-triplet structures 
were treated in a unified way. Further, in Ref|3J), which is referred to as (C), a phase 
diagram was given on the x-iV plane, where N represents a quark number. However, 
in the series of previous papers (A)~(C), a color neutral quark-triplet state was only 
realized ELS cL color-singlet state. 

In the first paper of the present series, namely in Ref|5]), which is hereafter 
referred to as (I) , the exact eigenstates are constructed so as to satisfy a certain con- 
dition which gives a color-singlet state in average. As a result, the color-singlet state 
is obtained in the color-symmetric form. In this paper, which is the second paper 
of the present series, the ground-state energy is reinvestigated under the condition 
giving a color-singlet state. In each region of the force strength X: the character of 
ground state is investigated, in which the quark-pairing state or quark-triplet state 
as a nucleon may be realized with a distinct value of a certain variable which is 
regarded as the order parameter of a phase transition. It is shown that there are 
two transition regions with different values of the order parameter. 

This paper is organized as follows: In the next section, the basic scheme for 
searching the ground-state energy is given in the modified Bonn quark model. In 
§3, the condition for the ground state is investigated and in §4, the ground state is 
determined and the ground-state energy is derived in each area of the force strength 
X- In §5, numerical analysis is given. The last section is devoted to concluding 
remarks. 



Our main concern in the present paper (II) is to search the minimum values of 
energies corresponding to the ground states. For the convenience of the discussion, 
we extract the relations (1-6-17) and (1-6-19): 



§2. Scheme for searching the ground-state energies 





El m) (N°,n ) 



\r = n ) 




\N°Un° + 6 - N°) + -N°(N° - 6) for the area A s , (2-2b) 
4 6 
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Here, the areas A\ and A s are shown in Fig[TJ Later, we will treat A\ by decom- 
posing it into A\ x and A\ 2 which indicate the areas related to < N° < 3i?° and 
3J?° < iV° < 4J?°, respectively. The relation ([22]) shows us that ^ (m) (iV , re ) and 

, re ) are quadratic in n° for a given value of N°. In each area of A[ and 
A s , regarding n° as continuous variable, we first calculate the minimum values of 
El m) (N°,n°) and Ei m) (N°,n°), respectively, for a given value of N° and we compare 
with each other the two minima so obtained. The above is our scheme for searching 
the ground-state energies. 

First, we notice that the term F = (1/2) • (1 + 2%)F(N , n°; 2r) in the expression 
(|2-ip is negative definite, if \ obeys 

1 + 2 X <0, i.e., X<~\- (2-3) 

Therefore, in this case, the maximum value of F(N , n°; 2r) determines the minimum 
value of energy, which appears at the point 2r = 0, that is, 2s 7^ or 21 7^ 0. This 
case is contradictory to the requirement that F(N° ,re°;2r) should be as small as 
possible. Next, we consider the case 

1 + 2 X = 0, i.e., X = ~\- (2-4) 

In this case, independently of the magnitude of F(N°, n°; 2r), the term F vanishes. 
This indicates that the energies in the case 2r < 2r m are the same as that in the 
case 2r = 2r m , in other word, the case 2r = 2r m is not toward the smaller direction 
in energy. For the above reason, we will be concerned only with the case 

1 + 2 X >0, i.e., X>~\- (2-5) 

Certainly, the procedure for minimizing the energy in the case (|2-5p automatically 
leads to the condition required to F(N°, n°; 2r). The above tells that, in the case 
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X < —1/2, the present model loses its meaning for the many-quark model. This 
point was already stressed in (A) qualitatively. 

In (A) and (C), we showed the minimum value of the energy as a function of 
N° and the value of nP which minimizes the energy. In (A) and (C), we treated the 
case no = 0, but essentially the same is valid for uq ^ 0. Also, we mentioned that n° 
can be regarded as the order parameter of the phase transition between the quark- 
pairs (n° = 0) and the quark-triplets (n° = A°/3). The results obtained in this 
investigation seem to be quite natural and acceptable. However, this investigation 
was not based on the "color-singlet" states minimizing F(N°,n°;2r). Further, the 
results were provided only for the area A\. Therefore, we must reexamine the results 
presented in (A) and (C). 

Following the procedure already mentioned, we are able to obtain the result in 
each area. In the area Ai, we obtain the following: 

(i) n° = for < N° < 3f2° , (2-6a) 
. . 1 12° + 6 1 
(Zi;2) -6'7^ <X --6' 

.0_ n ^ ,rO/6(2^-3-6 X ) 



(i) n u = for < iV u < 



5 + 6 X 



(ii) n° = £ for <KM°-3-6x) < < 3^0 (m) 
3 5 + 6x 

, . 1 2J?°-3 

Gi;3) -77<x< 



6 ^ ~ 6 

2J7° 

i) n° = for < N° < 



(ii) n' 



l + 2x 

_ (JV° - 2^2° + 3) + 2 X (A° + 3) 
1 + 6 X 



for J^_3 <A rO< 3 ^o_9_ 

1 + 2 X ~ ~ 2 A ' 

/V^ Q 

(in) n ° = for 317° — - — 9% < A° < 3/2° , (2-6c) 

3 2 

/i ^ 212° -3 

(ii;4) — - — <x<+°°, 



N° 

(i) n ° = — for < A° < 312° . (2-6d) 
3 

The results (|2-6|) coincide with those shown in the relation (C-4-1). In the case (|2-6b|) . 
we find that, at the point A = 6(212° — 3 — 6x)/(5 + 6%), a phase transition occurs. 
The order parameter n° changes from n° = to n° = N /3. In the area Aj 2 , we 
obtain the following results: 

1 1 2^2° + 3 

(/ 2 ;1) -- < x <-_. 4^0 + 3 ' 
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(i) n° = for 312° < N° < 4f2° , (2-7a) 
1 2Q° + 3 1 
(/2;2) "2"4?^T3 <X -"6 ' 

(i) n0 = for 3^<jV0< 2 ^°- 3 ) + 6 ^ - 1 ^, 
U " 3 + 10x 

(ii) n° = 4^°-iV 

for 2[(4^-3) + 6(2^-l) X ] 0<4J?0 

3 + 10x " 1 ; 

(fa; 3) -^<x<+oo, 

(i) n° = 4tt° - N° for 3/2° < iV° < 4J?° . (2-7c) 

In the case (i2^7b|) . the ph ase transition occurs at the point N° = 2[(4^° - 3) + 
6(2,1?° — l)x]/(3 + 10x)- The parameter n° changes from n° = to n° = 4.1?° — N°. 
In the area A s , we have the following: 

, s _ 1 1 412° + 9 

lS ' j 2 <X " 6 '2^0 + 3' 

(i) n° = ^ for 0<jV o < 3(2^ + 3 + 6x) ; 

3 1 - 6% 

(ii) n° = iV° for 3(2^° + 3 + 6 X ) < < 2 ^o (2 . 8a) 

1 - 6x 

. 1 4^2° + 9 

(S;2) -6-27^T3 <X<+ °°' 
AT 

(i) n ° = — for < iV° < 2Q° , (2-8b) 
3 

(s;3) - - < x < +oo , 
N° 

(i) n ° = — for 2f2° < N° < 3f2° . (2-8c) 
3 

In the case ()2-8ap . a phase transition occurs at the point N° = 3(2i?° + 3 + 6%) /(l — 
6%). The parameter n° changes from n° = A^°/3 to n° = iV°. 

As was shown in (A), the present model can be formulated not only from the 
side N = but also from the side N = 61?. We called this form the hole picture, 
which is obtained by replacing iV° with (6,1?° — iV°) in the relations appearing in the 
form from the side N = 0. Therefore, it may be necessary to compare the results 
for A\ 2 with those for Ay obtained by replacing iV° with (61?° — N°) in A\ x . In the 
area Ay , we obtain the following: 

, 1 1 [2° + 6 

-2 <x -- D -7^' 

(i) n' Q = for 3f2° < N° < 61?° , (2-9a) 

. i n° + 6 i 
(«i;2) -6-7^<x<-g, 
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(i) n ' Q = 2Q -— for 3^ u < iV u < r 

3 5 + 6x 

= ta 1803° + 1)(1 + 2.Q < ^ < 6fl „ 
5 + 6x 

„/ \ 1 2J?°-3 
« ;3 » 

TV" 9 

(i) no = 2^° for 312° < iV° < 3^° + - + 9 X , 

3 2 

, _ (4tf° - iV° + 3) + 2x(6^° - iV° + 3) 

for 3^ + ^ + 9 X <jV°< 4 ( 1 + 3 ^° + 3 1 
2 A ~ ~ 1 + 2x 

(iii) n' = for 4(1 1 + + 3 ^° + 3 < iV < 6^2° , (2-9c) 

(q;4 — t — <X<+oo, 
o 

(i) n' = 2O - — for 3O < N° < 6J?° . (2-9d) 

§3. Condition for the ground-state energy 

In §2, we listed the condition for minimizing the energy in each area. In order 
to search the ground-state energy, we must compare the minimum energy value in 
each area with the other areas. First, we discuss the relation between and A s . 
As can be seen in the area A s , the relation n° = N° appears in the case (s; 1) (ii): 

2 A - 6 2f2° + 3 1 - 6x ~ 

In any other case, we have n° = N° /3. On the other hand, we notice the case (Ji; 1). 
Since -(1/6) • (4f2° + 9)/(2fl° + 3) < -(1/6) • (ft + 6)/(ft° + 2) and 212° < 3J?°, 
£^ m) (iV , ra° = JV°) should be compared with E\ m) {N°,n° = 0). For the comparison, 
we have the relation 

AE = #M (N°,n° = N°) - E\ m) (N°, n° = 0) 

= I ( 2^o _ jyO) + 1 (1 + 2x)(iV o + 2) (3 . 2) 



The relation (|3-2|) indicates that, in the case (|3- 1|) . Z\i? is positive. In any other case 
in A s , we have n° = N°/3 which belongs to the border of A s and A\ and we notice 
the relation 

E^ m \N°,n = N°/3) = E^ ] (N°,n = N°/3) . (3-3) 

Therefore, the comparisons of Ei m) (N°, n° = N°/3) with E[ rn) (N° , n°) is reduced to 
those of E^iN ,^ = N°/3) with E^(N°,n°): The results are reduced to those 
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shown in the relations (|2-6ap ^ (|2-6dp . There does not exist any contribution of the 
area A s to the ground-state energy. 

Next, we consider the overlap of the areas A[ 2 and Ay . The overlap in this case 
consists of the regions composed of the following combinations: 

d ; n° = , n 0/ = , (3-4a) 
C 2 ; n° = AQ° — N° , n 0/ = , (3-4b) 

C 3 ; n° = , n ' = 2f2° - — , (3-4c) 
C 4 ; n° = 4n° - N° , n 0/ = 2^° - — , (3-4d) 



C 5 ; B o = 4flO-^, w „ = (4tf , -iV + 3) + 2x(6fl°-iV0 + 3) j 

1 + 6x 

For the combinations Ci ~ C4, we can prove that AE{= E[ 2 — Ey) are positive: 

AE >0 for Ci ~ C 4 • (3-6) 
For this proof, the following forms are useful: 

AE7 = (l + 2x)(ft° + l)X for Ci , (3- 7a) 

^ = 1(3 + i 0x )X 2 + 1(1 + 2 X )(^° + 5)X 

-1 [(1 + 6 X )^° + 2(1 + 2 X )] fl° for C 2 , (3-7b) 

AE = 4(5 + 6 X )X 2 + 1(1 + 2 X )(^° + 1)X 
3d 2 

+1 [(l + 6x)^° + 6(l + 2 X )] f2° for C* 3 , (3-7c) 

AE = ^(1 + 3 X )X 2 + 2(1 + 2 X )X for C7 4 . (3-7d) 



Here, X denotes 



X = N° - 3f2° . (3-J 



For the combination C5, a rather lengthy consideration may be necessary. Our 
task is to investigate the overlap of the cases (I2; 3(i)) and (l[; 3(ii)) shown in the 
relations (|2-7c[) and (|2-9cp . respectively. In this case, the overlap range of X is, at 
most, 

1 2Q° - 3 , , 

-e <x< -5-- (3 ' 9 > 

Further, we have the relation 

3f?° + \(1 + 2 X ) < N° < min + 3, 412°) . (3-10) 
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(III)- 









(II) 
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Fig. 2. The behavior of <?s(x)> 9d(x) an d <?o(x) are illustrated. In the area (I) and (II), AE < 0. 
The area (III) is divided into two areas. 



Of course, the relation (|3-10p gives us 
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3^° + -(l + 2x) <4tt° 



The relation (|3Tip is rewritten as 

n°>g B ( X ), g B (x) = l(l + 2x) , i.e. 



77 < X < 

D 



2f2° -9 
~18 



The relation (|3-10p leads to the following relations: 

(tt° + 3) + 6(^° + l)x 



9b(x) < X < 



g B ( X ) < X < [2° . 
The relations (|3- 13|) and (|3-14p are valid, respectively in the ranges 

1 3 
~6 <X< ~2(2J?° + 3) ' 



2(2/2° + 3) 
The energy difference AE is given as 



< X < 



2f2°-9 
~18 



(3-11) 

(3-12) 

(3-13) 
(3-14) 

(3-15) 
(3-16) 



AE = 4x(l + 3x) x2 + (l + 2 X )(l + 12x) jy 



+ 



1 + 6x 1 + 6x 

6x(l + 2x)^ 0+ 9 (l + 2x) 2 



l + 6x 



4 l + 6x 



Since X > 0, all the terms appearing in the function f{X) are positive, if 
Therefore, we have 



(3-17) 
X > 0. 



if x > , AE > in the range X > 



(3-18) 
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Fig. 3. The behavior of gi(x) an< A 9z{x) are illustrated together with <?b(x)> S-d(x) an< A 9o(x)- The 
area (III) is divided into two areas, (III)' and (III)". In the area (III)', AE < and in (III)", 
the regions of AE < and AE > coexist. In the area (IV), AE < 0. In the area (V), AE < 
and in (V)', the regions of AE < and AE > coexist. Also, in the area (VI) and (VI)', the 
regions of AE < and AE > coexist. In the area (VII), AE > 0. The area (V)' is too small 
to show the size definitely. 



The above means that if x > 0, the energies in Ai 2 are larger than those of Ay for 
N° > 3i?°. From the above reason, we will investigate x < in the cases (|3- 13|) and 
(|3-14p . together with the relations ()3- 15j) and (|3T6p . Of course, the relation (J3T6P 
is changed to -3/(2(21?° + 3)) < \ < 0. Then, the relations (13T5]) and (^46]) are, 
respectively, rewritten as 

n° > g (x) , (3-19) 
Q° < g ( X ) , (3-20) 

9o(x) = - r —^. (-6<X<0j (3-21) 

First, we treat the case (|3-13j) combined with the relations (|3-12p and (j3- 19|) . 
The functions related to this case are gsix)-, 9o(x) an d 9d(x)- The functions gs (x) 
and go(x) are defined in the relations (|3-12p and (|3-2ip . respectively. The function 
9d(x) is defined through the discriminant of the quadratic function f(X) given in 
the relation ([gT7jl . D: 

D ._?!^a+W±M (fl ._ OT(x)) . (3 . 22) 
9D[X] = C + W-ax)' {m 

The behaviors of <7s(x), go{x) an d 9d(x) are illustrated in Fig(2j Since gB(x) > 
and <7o(x) > 0; < an d ^-E > are determined by f2° > <7d(x) and i?° < 9d{x)i 
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Fig. 4. The detail around the area (V)' in Fig. 3 is revealed. 



respectively. In the areas (I) and (II), AE < 0. In (III), D > 0, but, in the present 
frame, AE > or < cannot be determined. 

Next, we treat the case (|3-14|) combined with the relations (|3-12|) and (|3-2U|) . 
In addition to 0b(x)> 9o(x) an d gD(x), we introduce the functions gi(x) an d gi (x) 
which are defined as follows: 

fiSBtx)) = ^iVg^ {f] ° ~ 9lU)) = h ' (3 ' 24a) 

3l(x) = -^-(l + 2x) 2 (l + 18x) , (3-24b) 
8x 



l + 6x V 16x(l + 3x) 92(x)y l 7 1 ^ 

g 2 (x) = - g^!^^) (! + 18 X + \A + 216x 2 ) • (3-25b) 

Here, it should be noted that (J2° - 9(1 + 2x) 2 /( 16 x(i + 3x)) ■ (1/02 (x))) is positive. 
It may be clear from the definition of gi{x) for i = 1 and 2, that if 12° < gi(x), fi > 
and vice versa. 

Figure [3] shows the behaviors of <7i(x) an d <?2(x)> together with 5s(x) ; <?o(x) 
and 0£)(x) which are shown in Figj2j With the use of the relation between gt{x) 
and /j, we obtain the following feature for AE < and AE > 0: (i) In the area 
(IV), AE < 0. (ii) In the area (V), AE < and in (V)', the regions of AE < 
and AE > coexist, (iii) In the areas (VI) and (VI)', the regions of AE < and 
AE > coexist, (iv) In the area (VII), AE > 0. Detail discussion on the coexistence 
of AE < and AE > is done in next section. 
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Finally, we discuss the area (III) in Figj2j As is shown in Fig|3j the area (III) 
is divided into two areas, (III)' and (III)". In the area (III)', AE < and in (III)", 
AE < and AE > coexist. In next section, we will discuss the meaning of the 
linex = -(4-v / 13)/18. 

§4. Determination of the ground-state energies 

We will continue the discussion in §3. The present model contains two param- 
eters, £2° and x- The parameter i7° denotes the available degeneracy of the single- 
particle levels and determines the framework of the model. On the other hand, x 
determines the dynamics caused by the model. Therefore, it may be natural to inves- 
tigate the dynamics by changing x f° r a fixed value of From Figs J2] and [3l we can 
learn that if x changes from x = ~ 1/6 to x = for a given there appear various 



features. The case 9/2 < fl° < Q° c = (7 + v / 13)/3 + ^3(28921 + 4784 v / 13)/54(= 
10.42 • • • ) may be the simplest, but physically interesting. Other cases are trivial or 
too complicated. In this paper, we will investigate mainly this case. 

Originally, f2° is a positive integer and in the present case, J2° takes the values 

as 

[2° = 5, 6, 7, 8, 9 or 10. (4-1) 

In (A) and (C), we showed various numerical results in the case J?° = 6. Therefore, 
the comparison of the present results with them may be interesting. We can classify 
the range — 1/6 < x < into the following cases: 

(1) - \ < X < XoW ) m (I) , (4-2a) 

b 

(2) xo(^°) < X < in (IV) + (V) , (4-2b) 

(3) xi(^°)<X<X2(^°) in (VI) , (4-2c) 

(4) x 2 (^°) < X < in (VII) . (4-2d) 

Here, Xi(" ) (* = 0> 1> 2) denotes the inverse of f2° = Qi{x)- The function xo(^°) 
is simply expressed as 

The function xi(" ) is obtained in an extremely complicated form: 

Xi(n°) = U/-b + VR + \J-b-VR - 
Z° = 4f2° + 99 , R = a 3 + b 2 , 

a = — — (Z 02 - 4617) , b = — ^— (2Z 03 - 13851Z + 177147) . (4-3b) 
729 v ; 19683 v ) \ J 

The function X2(^°) is expressed in the form 

= 4 4 "° + 9 , . (4.3c) 

V ; 2 (2tt° + l)(2fl° + 9) + 2f2°V^ 02 + 281?° + 55 
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For the case £2° = 6, we have 

Xo (6) = -0.1 , xi (6) = -0.041099 , xa(6) = -0.032811 . (4-3d) 

With the aid of the properties of the quadratic function f(X) (in — oo < X < 
+oo), we have the following features: 

(1) ' AE < for , flW < X < (tf+3) + 6 (fl ° + l)x 

1 + 2x 

(2) ' AE < for ffs (x) < X < Q Q , (4-4b) 

(3) ' AE > for g B ( X ) <X< X L ( X ; , 

AE < for X L (x; f2°) < X < f2° , (A- Ac) 

(4) ' AE > for g B {x) < X < f2° . (4-4d) 

Of course, (1)' ~ (4)' correspond to (1) ~ (4), respectively. Here, Xi denotes the 
larger for two solutions of f(X) = and it is explicitly given as 

Here, f(X) is given in the relation (|3-17p . The equation f(X) = is also regarded 
as quadratic for \ an d the larger solution is obtained in the form 

1 AX + 9 

X = 



2 2(3^° + X{2X + 5)) + (AX + 9) + 2^(3^2° + X(2X + 5)) 2 - X 2 (AX + 9) 
= Xl(X;I2°). (4-5b) 

The first term in the relation (|4-5ap denotes the value of X which makes f(X) 
maximum and it is smaller than gsix) m the case 

4-\/l3 N 
X < ^— ' ( 6) 

Therefore, the smaller solution loses its meaning. 

Until the present, we have investigated the effect coming from the area A[ 2 in 
the range 3i?° < iV < Af2°. By replacing N with (<6fi — N), the above effect can be 
included in the range 2f2° < N° < 3f2°. Including other areas, we can arrange the 
results in the ground-states as follows: 

m -I -I ^° + 6 

lj 2 <X< 6'tt° + 2' 

(i) < N° < 3Q°, n° = , (4-7a) 

fn , 1 f2° + 6 1 

(2) -6-7^ <X< "6' 

. 6(212° - 3- 6x) n 

i < iV° < 1 — ^ , n° = , 

5 + 6x 
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(ii) 6(2 ^; 6X) <^<3^, nO = ^, (4.7b) 
(3) -\< X <Xo{^), 

(i) 0<iY° <T ^--3, n° = 0, 

2S ' >] -3<N <3n -l(l + 2 X ), «°-«°t 



(4) 



Q A7^0 

ii) 312° - -(l + 2x) < iV° < 3^° , n° = — , (4-7c) 
2 3 

Xo(^°)<X<Xi(^°) , 
(i) 0<iY0 <T ^--3, n° = 0, 
2^° 



( 



1 + 2 X 2 



n = n 



3 < N° < 2f2° , n° = n°* , 

1 + 2x 

ii) 2(2° < N° < 3n° - H(i + 2 X ) , n° = n ot 



IV 



3J?° - -(1 + 2 X ) < N° < 3Q Q , n° = , (4-7d) 

(5) xi(^°)<X<X2(^°) , 

(i) 0<N»<-^--3, n° = 0, 

_ 3 <iV <2^ , n° = n *, 



ii) 212° <iV° <3^°-X L (x;^°) , n° = n°t 
iv) 312° - X L { X ; f2°) <N °< 3^-^(1 + 2 X ) , 



n° = n°* 



Q AT® 

(v) 3^° - -(1 + 2 X ) < N° < 3f2° , n° = — , (4-7e) 

(6) X2(^°) < X < , 

(i) 0<N°<-^--3, n° = 0, 

~ —0 ,,n0 9, . o „0* 



(ii) — - 3 < N° < 3f2° - -(1 + 2 X ) , 



n = n 



1 + 2 X ' 2 

q jv° 
(iii) 312° - |(1 + 2 X ) < iV° < 3^° , n° = — , (4-7f) 



_ 2tt° -3 

(7) < X < +oo , 

A/" 

(i) < N° < 3^° , n° = — . (4-7g) 
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Here, n ^ and n°* denote 

B t=jV-20», „0. = (iV°-2^ + 3) + 2^° + 3) 

1 + 6 X 1 ' 

If n ^ is replaced with n *, the result (|4-7p reduces to the result (C-4-1). With 
the use of the expression (|2-2ap . the ground-state energies are calculated except the 
case n° = n9\ in which the energies are calculated by the relation (|2-2ap under the 
replacement of N° with (6J?° — iV ). 

§5. Numerical analysis 

In this section, we give numerical results in the region A^. The quantity n , 
which is regarded as an order parameter of a phase transition between the quark- 
triplets and the quark-pairs, is given in Eq. (|4-7p in the various areas. The ground- 
state energies are also calculated numerically by using Eq. (|2-2a|) . In this section, we 
fix the parameter J2° = 6. Thus, we show the behaviors of the ground-state energy 
and the order parameter in the region < iV < 317° (= 18) with various force 




Fig. 5. The behaviors of energy (left panel) and the order parameter n° (right panel) are shown as 
functions of N° in the case (1) in Eq. (l4-7a|l with \ = —1/3. The model parameter f2° is taken 
as 6. 




Fig. 6. The behaviors of energy (left panel) and the order parameter n° (right panel) are shown as 
functions of TV in the case (2) in Eq. (l4-7b|) with \ = —1/5. The model parameter (2° is taken 
as 6. 
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5 10 15 20 

Fig. 7. The behaviors of energy (left panel) and the order parameter n° (right panel) are shown as 
functions of N° in the case (3) in Eq.(|4j7c} with x = —1/8. The model parameter Q° is taken 
as 6. 




Fig. 8. The behaviors of energy (left panel) and the order parameter n° (right panel) are shown as 
functions of iV in the case (4) in Eq. (l4-7d|) with x = —1/16. The model parameter J7° is taken 
as 6. 



strength of the particle-hole interaction, x- 

In Figj5j the ground-state energy on the left panel and the order parameter on 
the right panel are, respectively, shown as a function of N° with \ = —1/3, namely 
the results are based on Eq. (|4-7ap . In this case, the order parameter n° is identical 
to 0. Thus, the quark-pair state is realized. Almost the same behavior is seen in 
Fig.2 in (A). 

In Figj6j the ground-state energy and the order parameter are shown in the 
case x = ~ 1/5, namely the results are based on Eq. (|4-7bp . In this case, the order 
parameter n° changes from to N°/3. Thus, as the particle number N° increases, 
the quark-pair state changes to the quark-triplet state directly. This behavior is seen 
in Fig. 3 in (A). 

However, in FigJTl another behavior of the order parameter n° is seen. The 
ground-state energy and the order parameter are shown in the case x = —1/8 based 
on Eq. (|4-7c|) . In this case, the order parameter n° changes from to N°/3. However, 
in the region 13 < iV° < 14.625 (if iV° is integer, 13 < iV < 14), the transition region 
from quark-pair state to color-singlet quark-triplet state is open with n° = The 
change of order parameter is not continuous, while the order parameter is changed 
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Fig. 10. The behaviors of energy (left panel) and the order parameter n° (right panel) are shown 
as functions of iV° in the case (5) in Eq. l|4-7ff ) with \ = 1. The model parameter O is taken as 
6. 



continuously shown in Fig. 4 in (A). For the point between n° = or n° = iV /3 and 
n° = n ^ , the ground-state energy is not continuously changed with respect to the 
change of N°. The dotted curve in the left panel represents the ground state energy 
with n° = re°t . However, as the particle number iV° increases, the quark-pair state 
changes to the quark-triplet state through the transition region. 

In FigJSJ the ground-state energy and the order parameter are shown in the case 
X = —1/16 based on Eq, (|4-7dj) . In this case, the order parameter n° changes from 
to N°/3 through n° = n°* and n° = n ^ . For the point between n° = n°* or 
n° = N°/3 and n° = nP\ the ground-state energy is not continuously changed with 
respect to the change of iV . The dotted curve in the left panel represents the ground 
state energy with n° = n ^ . In this case, also, as the particle number iV° increases, 
the quark-pair state changes to the quark-triplet state through two transition regions 
characterized by n°* and n ^ . 

Also, in FigH] in the case % = —1/27 based on Eq. (|4-7e|) , the order parameter 
n° changes from to iV°/3 through n° = n *, n° = n ^ and n° = n°*. For the 
point between n° = n°* and n° = n°\ the ground-state energy is not continuously 
changed with respect to the change of iV similar to the case \ = —1/8. Here, the 
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Fig. 11. The behaviors of energy (left panel) and the order parameter n° (right panel) are shown 
as functions of N° in the case (5) in Eq.( |4-7g] l with x — 3. The model parameter Q° is taken 
as 6. 



dotted curve in the left panel represents the ground state energy with n° = n°t. In 
this case, also, as the particle number N° increases, the quark-pair state changes to 
the quark-triplet state through three transition regions. 

In Fig JIOl the ground-state energy and the order parameter are shown in the case 
X = 1, namely the results are based on Eq. (|4-7fj) . In this case, the order parameter 
n° changes from to iV /3 through the transition region with n° = n°*. Thus, as 
the particle number iV° increases, the quark-pair state changes to the quark-triplet 
state through the transition region. This behavior corresponds to Fig. 5 in (A). 

Finally, in Fig(TTJ the ground-state energy and the order parameter are shown 
in the case x = 3, namely the results are based on Eq.( |4-7g' ). In this case, the order 
parameter n° is identical to N°/3 in all the regions. Thus, the color-singlet quark- 
triplet state is realized. This behavior of the order parameter corresponds to that 
shown in Fig. 6 in (A). 



§6. Concluding remarks 



In this paper, the ground-state energies were calculated in the modified Bonn 
quark model. The color-singlet condition was imposed, which was developed in 
the previous paperA' It was indicated that the modified Bonn quark model has a 
meaning in the regions x > — 1/2, where x represents the strength of the particle- 
hole-type interaction. The ground-state energies were determined in each area of x 
and the ground states were also determined. Namely, in a certain parameter region, 
the color-singlet quark-triplet state is the ground state, and in another region, the 
quark-color-pairing state is the ground state. Further, it was shown that there were 
two transition regions between the quark-pairing and the quark-triplet state, which 
were distinguished by the value of the order parameter. 

Finally, we give a comment on the phase transition induced by the present model. 
For this aim, we show the phase diagram obtained by the relation (|4-7p . It is drawn 
for the case i?° = 6 in Fig]121 Since n° plays a role of the order parameter, the present 
model induces four phases and the areas specified by n° = N°/3 and n° = are in the 
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quark-triplet and the quark-pair phase, respectively. In the original Bonn model, the 
case x = was adopted. The region N° ~ or 612° and the region N° ~ 3J?° are in 
the quark-pair and the quark-triplet phase, respectively. Further, we must note that 
in the former and the latter regions the quark system has a low and a high density, 
respectively. We observe, however, that the situation predicted by the present model 
is in contradiction with the common understanding p» according to which color- 
superconductivity is expected to provide a reliable description of quark matter at 
high densities. In the case where x is different from x = 0, the situation does not 
change, if we treat x as a constant in the whole region < N° < 612°. Although the 
original and the modified Bonn quark model present us various features of many- 
quark system, it may be not permitted for these models to avoid the investigation 
of the above-mentioned problem. In (III)]-" we will discuss this problem. 
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